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ABSTRACT:  A  class  of  random  transient  signals  has  been  defined  as 
the  product  of  a  deterministic  envelope  waveform  of  finite  integral 
square  and  a  continuous  random  process  with  a  well-defined  power 
spectrum  and  autocorrelation  function*  The  time  average  autocorrela¬ 
tion  function  and  energy  density  spectrum  of  the  resulting  waveform 
have  been  found  to  be  random  variables  at  every  value  of  their 
arguments.  The  means  and  variances  of  these  random  variables  are 
derived  as  functions  of  the  characteristics  of  the  envelope  and 
original  noise  process.  The  average  autocorrelation  function  is 
found  to  be  the  product  of  the  autocorrelation  functions  of  envelope 
and  noise,  and  the  average  spectrum  is  given  by  the  convolution  of 
the  energy  spectrum  of  the  envelope  function  and  the  power  spectrum 
of  the  noise.  Examples  of  the  mean  and  variance  calculations  are 
presented  for  both  rectangular  and  decaying  exponential  pulses  of 
both  broad  and  narrow  band  noise.  Finally,  the  implications  of  these 
findings  for  measurement  programs  and  monopulse  signal  processing  are 
discussed. 
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A  class  of  random  signals  has  been  modeled  as  the  product  of  a  trans¬ 
ient,  deterministic  envelope  waveform  and  a  well-behaved  continuing 
random  process.  The  properties  of  the  energy  density  spectrum  and 
autocorrelation  function  of  such  signals  are  studied  and  the  results 
related  to  current  problems  in  signal  processing  and  monopulse  detec¬ 
tion  systems.  The  work  on  this  project  was  funded  under  Task 
ASW2-21-000-W270-70-00.  The  report  will  be  of  interest  to  those 
concerned  with  statistical  communication  and  detection  theory,  active 
sonar  systems,  signal  processing,  and  noise  immunity  studies. 
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Chapter  I 
INTRODUCTION 


In  the  study  of  a  large  class  of  communication  and  signal  detec¬ 
tion  systems,  one  is  often  faced  with  the  analysis  of  the  effects  of 
interfering  noise  of  a  transient,  non-continuing  nature.  Examples 
of  such  noise  phenomena  include  the  reverberation  background  of  a 
sonar  signal  and  impulsive  interference  of  the  type  seen  on  telephone 
lines  and  atmospheric  radio  links.  When  formulating  a  system  for  the 
detection  of  wanted  signals  in  such  a  background,  it  is  often  neces- 
ry  to  know  in  some  detail  the  frequency  distribution  of  energy  in 
,he  interference  or  how  such  a  noise  coraoonent  behaves  under 
correlation  processing.  The  intent  of  this  report  is  to  detail  an 
investigation  of  certain  statistical  properties  of  a  class  of  noise 
bursts  suggested  by  the  above.  It  is  hoped  that  the  results  gained 
here  place  in  somewhat  better  perspective  the  problems  faced  in  the 
analysis  and  synthesis  of  processing  systems  working  against  non¬ 
stationary  backgrounds. 

The  noise  signals  to  be  treated  here  are  neither  continuing 
stochastic  processes  in  the  usual  sense  nor  deterministic  transients 
amenable  to  immediate  treatment  by  the  Fourier  integral.  They  share 
the  properties  of  both  broad  classes  but  lack  the  mathematical 
convenience  that  arises  from  the  usual  assumptions.  (It  is  suggested 
that  these  signals,  bearing  many  properties  of  both  transients  and 
random  signals,  be  known  as  "random  transients".)  Since  a  noise  burst 
is  defined  only  for  a  given  epoch,  its  statistical  properties  are 
tied  to  a  given  instant  of  time,  and  stationarity  disappears.  Since 
ensemble  averages  no  longer  equal  time  averages,  ergodicity  soon 
evaporates  also.  On  the  other  hand,  such  a  burst  does  not  have  a 
deterministic  Fourier  transform,  and  Fourier  integral  analysis  must 
be  approached  with  great  care.  Even  so,  by  carefully  defining  terms 
and  remaining  reasonably  aware  of  the  necessity  of  continually  relat¬ 
ing  the  mathematics  to  the  physical  situation,  it  is  possible  to  achieve 
a  consistent  and  useful  interpretation  of  impulsive  noise  phenomena. 

Eventually,  it  emerges  that  the  autocorrelation  functions  and 
spectra  of  such  noise  signals  no  longer  possess  deterministic  values 
at  every  point,  but  rather  become  random  variables  with  calculable 
means  and  variances.  Fortunately,  it  is  possible  to  show  that  the 
means  are  given  by  quasi-intuitive  expressions  similar  to  those 
developed  in  traditional  transient  or  random  theory.  The  variances, 
in  turn,  provide  an  indication'  of  the  uncertainty  of  the  spectra  and 
correlation  functions  for  a  given  value  of  their  arguments.  The 
analysis  thus  lends  a  good  deal  of  insight  to  noise  measurement 
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programs  and  also  to  the  choice  of  a  processing  scheme  that  allows 
sufficient  latitude  to  encompass  the  great  majority  of  interfering 
background  noise  that  may  arise. 
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Chapter  II 

A  MODEL  FOR  THE  GENERATION  OF  RANDOM  TRANSIENTS 


The  model  used  to  generate  the  transient  noise  signals  to  be< 
treated  herein  is  portrayed  in  Figure  1.  The  random  transient  is 
taken  as  the  output  of  a  multiplier  whose  inputs  are  a  zero-mean 
Gaussian  random  process  n(t)  and  an  "envelope  waveform"  e(t),  both 
of  which  are  real.  The  former  is>  assumed  to  be  stationary  and  ergodic, 
thus  possessing  a  well-defined  power  spectrum  and  autocorrelation 
function  as  described  in  reference  (a).  The  "envelope  waveform"  is 
held  to  be  a  deterministic  transient  of  finite  integral  square  which 
is  zero  for  t<0.  The  output  of  the  multiplier  is  the  product  of 
these  functions  and  evidently  equals  zero  when  e(t)  ■  0. 

Some  justification  of  this  model  is  provided  by  noting,  that  many 
of  the  processes  which  yield  random-transient-like  signals  can  be 
approached  from  a  theoretical  basis  which  yields  a  prediction  of  some 
quantity  which  can  loosely  be  described  as  the  "average  level"  as  a 
function  of  time.  In  sonar  applications  for  example,  it  is  possible 
to  derive  theoretical  expressions  for  the  acoustic  power  returned  as 
a  function  of  range  from  either  volume  or  boundary  reverberation. 
Similarly,  in  the  study  of  transients  caused  by  impulsive  phenomena 
such  as  chemical  explosions  or  spark  gaps,  a  theoretical  treatment 
may  well  provide  an  expression  for  some  kind  of  envelope  within  which 
the  detailed  structure  of  the  transient  is  more  or  less  random.  It 
is  certainly  visionary  to  claim  that  such  a  multiplicative  envelope 
function  can  be  rigorously  defined  for  the  physical  processes  of  the 
real  world..  The  point  is,  however,  that  although  the  detailed 
structure  of  a  particular  transient  may  be  vastly  different  from  every 
other,  such  a  .quantity  as  an  average  "level"  or  perhaps  some,  analogous 
statistical  measure  may  well  display  a  uniformity  from  sample  to 
sample  that  can  be  described  with  the  artifice  of  a  multiplicative 
envelope. 

Here,  this  envelope  waveform  has  been  considered  a  deterministic 
transient,  probably  the  simplest  assumption  that  might  haye  been  made. 
It  is  felt  that  such  an  approach  is  consistent  with  the  class  of 
situations  described  immediately  above,  but  certainly  a  more  elabo¬ 
rately  structured  model  could  be  envisioned  to  accomodate  a  larger 
class  of  natural  phenomena.  As  a  first  step  upward  from  the  present 
assumption,  one  might  randomize  the  envelope  by  considering  such  para¬ 
meters  as  length  and  amplitude  to  be  random  variables  with  appropriate 
distributions.  Another  possibility  is  to  model  bhe  envelope* waveform 
itself  as  a  segment  from  a  relatively  low  frequency  process  so  that 
both  functions  entering  the  multiplier  are  completely  random.  These 
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FIGURE  1.  IDEALIZED  MODEL  Or  WAVEFORM  GENERATION. 
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approaches  are  not  developed  here  but  remain  g^iSiSStlSatSn”!? 
for  future  study.  At  least  at  the  outset  of  this  investigation 

was  felt  that  the  simplest  model  was  adequate  to  deal  with  atten- 
phenomena  of  immediate  interest.  We  shall  now  concentrate 
tion  on  the  model  of  Figure  1. 

The  term  "envelope  waveform"  is  actually  something  of 

slowly  varying  envelope  of  a  more  rapid  oscillation.  Consider  the 
output  waveform: 


s(t)  *  e(t)  n(t) 


where  at  every  instant  of  time  t,  n(t)  is  a  Gaussian  random  var¬ 
iable  with  mean  zero  and  variance  rj^  •  At  this  same  instant  of  time, 

N 

s(t)  is  thus  also  a  Gaussian  random  variable  with  zero  mean  and  with 
variance.  .  - . 

02  (t)  =  e2Ct)a2  (2) 

s  N 

Alternatively, 

as(t)  =  |e(t)|aN  j  (3) 

and  it  is  seen  immediately  that  the  most  meaningful  interpretation  of 
e(t)  is  as  the  time  varying  factor  whose  magnitude  relates  the  stand-  ' 
ard  deviations  of  the  input  and  output.  It  is  also  apparent  now  that 
s(t)  can  in  no  way  be  considered  a  stationery  random  variable  since 
indeed  its  variance  is  a  function  of  time.  This  is  hardly  surprising 
since  ve  are  now  dealing  with  a  transient  signal  for  which  concepts 
of  stationarlty  are  irrelevant. 

Consider,  however,  an  infinite  ensemble  of  waveform  generators 
on  this  model,  each  with  its  separate  independent  random  process 
n(t),  but  with  Mentical  envelope  waveforms  e(t).  The  ensemble  of 
output  transients  generated  by  such  an  assemblage  can  be  expected  to 
display  a  large  measure  of  statistical  regularity,  and  it  is  with 
this  ensemble  of  all  possible  transients  with  identical  e(t)  that 
we  propose  to  work. 
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CHAPTER  III 

NOISE  BURST  AUTOCORRELATION  FUNCTIONS 


For  a  real  deterministic  transient  such  as  the  envelope  wave¬ 
form  e(t),  the  autocorrelation  function  is  generally  defined  as 


%eW  5  fe 

4o< 


(t)e(t+?)dt 


such  that 


%e(0)  « ye2(t)dt 


where  See(0)  is  known  as  the  energy  of  e(t)  in  the  sense  that  if  e(t) 
were  a  voltage  or  current  waveform,  this  expression  gives  the  total 
energy  dissipated  by  e(t)  in  a  pure  one  ohm  resistance  (reference  (a)). 

For  a  real  random  process  such  as  n(t),  extending  in  time  from 
-<*>to  +oo  ,  two  autocorrelation  functions  can  be  defined.  First 
the  ensemble  average  autocorrelation  function: 


Rn  (^1**2^  a  E  n(t1)n(t2)J 


where  E[  ]  denotes  the  statistical  expectation  taken  over  all 
members  of  the  ensemble.  Rn  is,  in  general,  a  function  of  t]_  and  t2» 
but  if  the  process  is  stationary,  it  becomes  a  function  only  of  V 
the  difference  between  t]_  and  t2: 


Rn(ti,t2)  =  Rn(r)  -  E  |n(t )n(t+f)J 


where  X  =  t^-t2  and  t^  is  arbitrary 
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which  is  the  average  power  of  n(t)  in  the  sense  that  this  expression 
gives  the  average  power “dissipated  in  a  one  ohm  resistor  by  n  (t). 

If  the  process  is  ergodic,  time  and  ensemble  averaging  are 
equivalent,  and  in  particular, 

Rn(r)  =  fadr)  (10) 

With  this  background,  we  can  proceed  to  a  meaningful  definition 
and  evaluation  of  the  autocorrelation  function  of  the  signal  s(t) 
defined  in  equation  (1)  and  Figure  1.  Consider  first  the  attempt  to 
form  an  ensemble  average  autocorrelation  function  for  s(t)  as  in 
equation  (6): 

Rs(ti,t2)  =  E  [sUjXtaQ 

=  Efe(t1)e(t2)  n(ti)n(t2)J  (11) 

Since  e(t)  is  deterministic,  we  may  write 

Mti,t2)  =  e(ti)e(t2)E  [n(ti)n(t2j]  (12) 

Because  n(t)  has  been  assumed  to  be  ergodic,  this  becomes 

R5(ti,t2)  =  Rs(ti,T+ti)  =  e(ti)e(t!+r)  Ynn(f)  (13) 

r  =  t2-tx 

Since  Rs  is  a  function  of  both  ti  and  T  ,  s(t)  is  nonstationary,  and 
the  ensemble  average  autocorrelation  function  loses  much  of  its 
interest.  Such  is  generally  the  case  in  treating  a  transient  waveform. 


*  It  may  be  well  to  point  out  that  in  this  report,  the  letter  P (phi) 
will  be  used  to  denote  the  autocorrelation  functions  of  energy  signals 
as  in  equation  (4),  whereas  the  letter  f  (psi)  will  be  used  tor  the 
autocorrelation  functions  of  power  signals  as  in  equation  (8). 
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By  using,  however,  the  formula  of  equation  ( k )  which  is  appro¬ 
priate  for  traditional  transient  analysis,  we  can  write  that 


■  / 


sCtJaCt+Y)  dt 


e(t)e(t+r)n(t)n(t+T)  dt 


(14) 


where  the  lower  limit  of  the  integral  becomes  zero  since  e(t)  *  0 
for  t  <  O,  Now  since  n(t)  is  a  random  variable  for  all  t,  it  be¬ 
comes  apparent  that  the  integral  of  equation  (14),  if  it  exists,  is 
also  a  random  variable.  The  existence  of  an  integral  such  as  that  of 
equation  (l*f)  is  examined  in  reference  (b)  which  treats  expressions 
of  the  form 


y(T)  =  /  h(t,r)x(t)  dt 

A 

where  only  x(t)  is  a  random  variable.  Now,  if 

jf  E[|h(t,T)x(t-)j]  dt  s  £  |h(t ) |  E  Qx(t))]  dt 


(15) 


<  CMC* 


(16) 


then  y(T)  exists  for  all  sample  functions  of  x(t)  except  for  a  set 
of  probability  zero,  and  furthermore 


i[y(r)j  *  J  h(t,T)  E[x(t)]dt 


(17) 


Looking  at  equation  (l1*)  in  this  light  we  see  that 


Elje(t)e(t+r)n(t)n(t+r)T]  dt  »  /|e(t)e(t+T)|E  Qn(t)||  n(t+r)jjdt  (18) 
°  r  a 

Now  El  |n(t)|  |n(t+T)[J  can  be  interpreted  as  the  autocorrelation  func¬ 
tion  of  a  full-wave  rectified  version  of  n(t),  which  exists  if  n(t) 
is  a  sufficiently  well-behaved  stochastic  process  as  we  have  assumed. 
But  by  the  basic  property  of  the  autocorrelation  function  that  its 
maximum  occurs  at  the  origin  for  any  random  or  transient  function, 

E  [|n(t)|  |n(t+r)|J<E  |jn(t)|  (n(t)|J  ■  E[n2(t)J  =  Vnn(°)  (19) 
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Therefore, 


yje(t)e(t+^)l  IS  jjn(t)n(t+t9(|  dt  ^  ^n(0)  J (eCt)!  |e(t+r)|  dt 


(20) 


<%„(0)  %eW<  °° 


since  n(t)  is  assumed  to  convey  finite  power  and  since  e(t)  is  a 
transient  with  finite  energy. 

It  has  thus  been  shown  that  ^ss(r^  exists  for  all  members  of 
the  ensemble  except  for  a  set  of  probability  zero  (whatever  that  means), 
and  it  begins  to  make  sense  to  speak  about  the  autocorrelation  function 
of  s(t)  in  these  terms.  Since  ^53(7)400  for  all  r  , 


\sW  =  T /  s(t)s(t+r)dt  =  J^! i V. 


(21) 


s  0,  for  all  r 

and  this  form  of  the  time  average  autocorrelation  function,  appropri¬ 
ate  for  continuing  random  signals,  becomes  meaningless  here.  Thus 
s(t)  is  an  energy  signal  rather  than  a  power  signal  and  in  that  sense 
is  more  akin  to  a  transient  than  to  a  continuing  random  process. 

If  one  defines  the  autocorrelation  function  of  s(t)  in  the  form 
of  equation  (*»-),  the  result  is  a  random  variable  whose  mean  and  vari¬ 
ance  at  every  Y  must  be  related  to  the  correlation  functions  of  e(t) 
and  n(t).  Thus  if 


^  J s(t)s(t+*)dt  =  y*e(t)e(t+r)n(t)n(t+'r)dt 


(22) 


the  ensemble  average  becomes 


E  D°ss  (T)]  =  E  \f  e(t)e(t+r)n(t)n(t+r)dtj 


(23) 


and  using  equation  (17),  this  is 


/©o 

e(t)e(t+r)  E[n(t)n(t+r)]dt 


(2*0 
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Finally,  by  equations  (^),  (7),  and  (10). 


E  jj°  ss(r)J  =  'l'rm(r)  <>ee(r)  (25) 


It  becomes  evident,  then,  that  the  expected  value  of  the  autocorrela¬ 
tion  function  is,  for  every  T  ,  equal  to  the  product  of  the 
autocorrelation  functions  of  the  envelope  waveform  and  the  original 
noise  process,  and  is  therefore  similar  to  the  result  gained  in 
seeking  the  autocorrelation  function  of  the  product  of  two  independent 
stationary  random  processes,  whore  the  autocorrelation  function  of  the 
product  is  equal  to  the  product  of  the  autocorrelation  functions. 

Some  care  must  be  taken  in  the  interpretation  of  equation  (25). 
What  is  being  claimed  is  precisely  this:  Given  a  sample  function 
from  the  ensemble  generated  by  the  model  of  Figure  1,  we  calculate 
the  time  average  autocorrelation  function  as  in  equation  (22)  -  by 
multiplying  the  sample  function  by  a  shifted  replica  of  itself  and 
integrating  the  product  from  0  to  «  ,  If  this  operation  is  performed 
on  a  number  of  sample  functions  and  the  results  averaged  for  fixed  ?' , 
the  average  will  tend  to  the  expression  of  equation  (25)  as  the  number 
of  sample  functions  grows  large. 

If  r  3  0,  equation  (25)  provides  the 

9 

E[*ss<°>]  =  ^nnCO)  *ee<0)  (26) 


numerically  equal  to  the  product  of  the  energy  of  e(t)  and  the 
average  power  of  n(t)  (bui  bearing,  of  course,  the  dimensions  of 
energy).  This  is  to  say  that  if  we  form  the  integral 


v’ss(O) 


(27) 


for  a  number  of  sample  functions,  the  average  value  approaches 
E  [v>ss(o)  as  the  number  of  sample  functions  taken  grows  large. 

Next  to  be  considered  is  the  variance  of  the  random  variable 
tfggCr)  generated  by  the  correlation  process  of  equation  (22).  From 
elementary  probability  theory  it  is  known  that 


Var|>ss(rj]  =  e[*>sS(t)]  -  {e  [»ss(r|]  }  2 


(28) 


or  in  words,  that  the  variance  is  equal  to  the  mean  square  minus  the 
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square  of  the  mean.  The  mean  square  of  must  thus  now  he 

calculated.  The  square  of  v>ss(r)  is  given  by 


2  C** 

v>ss(r)  s  J  e(t)e(t+r)n(t)n(t+r)dt  J  e(x)e(x+T)n(x)n(x+?)dx  (29) 


OO 

e(t)e(x)e(t+r)e(x+r)n(t)n(x)n(t+r)n(x+r)dx  dt 


Taking  the  expectation  and  interchanging  with  the  integration  yields 


-OO^OO 

to)]  zJ  /e(t)e(x)e(t+r)e(x+Y)  E[n(t )n(x)n(t+r)n(x+rj]  dx  dt  (30) 


At  this  point,  for  the  first  time,  we  must  make  explicit  use  of 
the  fact  that  n(t)  is  a  Gaussian  random  process.  Indeed,  it  can  be 
shown  (c.f.,  reference  (c))  that  if  n(t)  is  a  zero  mean  Gaussian 
process,  then 


E[n(t1)n(t2)n(t3)n(tlf)3  =  E[n(ti)n(t2)]  fifn^Mt^i] 


(31) 


+  E[n(t1)n(t35  ECn(t2)n(tif)] 


+  E[n(t1)n(tif)]  E[n(t2)n(t3)] 

Now  using  ergodicity  and  equations  (6),  (7),  and  (10),  we  may  write 
that 


E jn(t)n(x)n(t+r)n(x+r)J  =  +%n(t-x)  +  ^ir.(t-x-r)  ^(t-x+r)  (32) 


Following  through, 
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CO  oo 

[%g(r)]  zj  j e(t)e(x) e.(t+r)  e(x+r)  fn£(r)dx  dt 

Q  O 

oo  oo 

l  /e< t)e(x)e(t+r)e(x+r)[  ^nn^t-x)  +  (jx  dt 


(33) 


The  first  integral  may  immediately  be  written  as 

tyrmW  I  e(t)e(t+r)dt  j  e(x)e(x+r)dx  - •^n^(S')<^|(r) 

jo  yQ 

but  by  equation  (25), 

%£<?)%§(•«  -  jBfVssWjP  i 

and  directly,  by  cancellation  in  equation'  (28), 

VarJS^-Cvj]  zf  f  e  ( t )  e  ( x )  e  ( t4\T)  e  (  x+T) 

*'o  *v> 


(3*0 


(35) 


|^nn(^  rX  J  +  ^nn Ct -x+T)  3dn ( t  -x-r|dt  dx 


(36) 


The  evaluation  of  this  horrendous  integral  can  ;be  somewhat  eased  by 
making  the  change  of  variables  u  *  t-x  and  t  .  ;t,  and  then  elimin¬ 
ating  x  from  the-  equation  (|Jacobiah|  a  1;  first  quadrant  area  of 
tx  plane  corresponds  to  the  upper  half  of  the  tu  nlane)0 


'Var  [VsS(r)]  .f  f  e.Ct.)  e(t-u)e(  t+T)  e  ( t+T-u) 

Jo 

r 

l^n^(t)  +  ( u+  r)  ( u -  r) J  du  d t 

Reversing  the  order  of  integration  then  yields 

ob 

Var|%s(r)]  =  nn^u)  +  %n(u+T)  ^(u-r)] 


(37) 


f 


(38) 


i 


e(t)e(t-u)e(t+T)e(t+7-u)dt  du 
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Consider  now  the  inner  integral, 

r- 

J  e(t)e(t+t')e(t-u)e(t+,r-u)dt  = 
0 

where 


which  can  be  expressed  as 
/pit  »r)p(t-u,f)dt 


p(t,T)  s  e(t)e(t+r), 


(39) 

(40) 


a  function  defined  as  e(t)  multiplied  by  a  shifted  replica  of  itself. 
p(tj't-)  possesses  an  autocorrelation  function  defined  as 


.  op 

^pp(u,t;  =  J  p(t,r)p(t-u,T)dt, 


(41) 


identical  with  the  inner  integral.  The  variance  may  now  be  written 

pp 


var  [>si<r)]  -fi  („,♦)  [y^(u)  ♦  du 

-oo  " 


(42) 


As  defined  above,  (u,  r)  is  merely  the  autocorrelation 
function  of  the  signal  generated  by  multiplying  e(t)  by  a  replica  of 
itself  shifted  T  seconds.  Several  examples  of  this  function  will  be 
computed  in  a  later  section,  but  some  general  properties  emerge  at 
once.  Since  v>pp(u,r)  is  a  bona  fide  autocorrelation  function  with 

u  as  the  independent  variable  and  T  as  a  parameter,  it  must  be  an 
even  function  of  u  and,  having  no  periodicity,  assume  its  maximum 
value  (for  given  ?  )  when  u  ■  0.  Also,  since  e(t)e(t+*')  is  identical 
with  e(t)e(t-r)  except  for  a  shift  in  the  time  origin,  ¥>pp(u,tr) 

is  an  even  function  of  the  parameter  T  and  for  a  given  u  obtains  its 
maximum  value  when  f  «  0.  Since  the  integrand  factor  in  square 
brackets  (equation  (42))  is  also  even,  the  variance  may  also  be 
written  as 


Var  C^ss^)] 


2 


/<*> 

*pp<utT)  [Vm(u)  +  Y'nnfa  +  'r)  ^(u  -r)J  du 


(43) 


The  expression  for  the  variance  derived  here  is  particularly 
easy  to  evaluate  on  a  digital  computer  using  numerical  integration 
techniques,  and  several  examples  of  this  calculation  will  be  pre¬ 
sented  later.  The  result  is  to  be  interpreted  as  supplying  tne 
variance  of  the  random  variable  ^ss(T^  as  a  function  of  r  . 

Since  by  definition 

2 

Var  [9ss<T>]  -  Ej[/ss(r)  -  B  0ss(t)]]  )  (MO 
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we  have  found  a  measure  of  how  tightly  the  distribution  of  %s(*y) 
.cleaves  to  the  mean,  and  hence  of  the  amount  of  variability  to" be 
expected  between  points  of  equal  on  the  empirical  autocorrelation 
functions  of  representative  samples.  Evidently,  when  Var 

becomes  small,  we  can  reasonably  expect  that  the  empirical  measure¬ 
ment  will  be  close  to  the  expected  value.  It  is  of  considerable 
interest  to  know  where  the  variance  is  a  maximum. 


Consider  then  the  expression  of  equation  (43); 

Joo 

<€pp(u,T-)^(u)dU 


CO 

K 


(45) 


(u,r)  %n(u +r)  %n(u-r)du 


* ■/  rnn 


Since  Hnn(u,0)i><fDD(u,/>'),  for  all  u,  then 


Tpp\u,wc.-*pp 

.OO  go 


f%p(uS 

•'o 


^Vnn^du  <, 


°)?*hmu)du 


(46) 


and  the  first  integral  is  obviously  a  maximum  for  7s  0.  It  is  also 
true  that 

f ^ p ( u ,?)  ( u+r) ^in(u-r) du  <f%p(u,0)  ^(u+r)  ^n(u-r)du  (47) 

*(0  o 

and  it  is  shown  using  the  Schwarz  inequality  in  Appendix  A  that  if 


C?pp(u,0)£0  for  all  u,  then, 

oo  r  00 

/^p(u,0)  ^(u+r)  V'nntu-rMu  <J<fpp(u,0)^(u)du  (48) 

From  equations  (47)  and  (48)  then 

Var  <  Var  K.C0)]  for  all  r  (49) 


subject  only  to  the  condition,  essentially,  that  the  envelope  wave¬ 
form  be  everywhere  positive.  The  maximum  variance  thus  occurs  at 
the  origin,  as  one  might  have  expected  intuitively,  and  its  value 
is  given  by 
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Max{var  |>SS(T)]}  =  Var[<f|s(0)] 


=  ^ J %jj(u)  f e2(t)e2(t-u)dt  du  (50) 

0  J0 

Since  ^ppUjT)  approaches  zero  for  all  u  as  T  grown  large,  it 
should  be  apparent  that 


Var  [%s(r)]— ►  0  as  Y— .  00  (51) 

One  more  general  result  will  now  be  discussed  having  to  do  with 
the  behavior  of  the  autocorrelation  variance  as  the  extension  in 
time  of  the  envelope  waveform  increases.  If  the  statistical  proper¬ 
ties  of  the  original  noise  process  remain  constant  while  the  envelope 
duration  grows  larger  (in  the  sense  that  a  rectangular  pulse  of  width 
T  or  an  exponential  pulse  of  time  constant  T  become  '‘longer"  as 
T  increases),  we  might  reasonably  expect  t  3  variance  of  the  auto¬ 
correlation  to  be  decreased  since  in  effect  che  integration-averaging 
process  of  equation  (2.?)  is  being  carried  out  over  a  longer  and 
ionger  period.  In  other  words,  by  taking  a  longer  and  longer  "piece" 
of  the  input  noise,  one  approaches  the  operation  of  equation  (8) 
which  yields  an  expression  of  zero  variance.  Actually,  since  the 

T1^\Tunii.Yar;*'ar‘<?0  r(;la^G<3  to  the  average  total  energy  (see  equation 
C50)),  the  variance  will  increase  (absolutely)  as  the  envelope  dura¬ 
tion  (and  hence  its  energy)  increases.  For  this  reason  it  is  more 
instructive  to  work  with  a  normalized  form  of  the  maximum  standard 
deviation.  We  form,  then,  the  ratio  of  the  standard  deviation  at 
T  ■  0  to  the  expected  value  of  Cfss(r)  at  r  =  0: 


H  =  Aar  [<fs5(0)] 
°  ~  E  [«f5s(0J] 


(52) 


This  ratio  expresses  the  standard  deviation  of  <f~c(0)  as  a  percent¬ 
age  of  the  mean  of  ^^(O)  and  thus  measures  the  extent  that  the 
distribution  of  clusters  about  the  mean  in  the  worst  case. 

From  equation  (50) , 


R 


0 


(53) 
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Consider  the  following  heuristic  argument*  As  the  width  of  the 
envelope  function  increases  (as  V* nn(u)  remains  the  same),  ^pp(u,0) 

becomes  "wider"  also,  and  in  most  cases  can  be  considered  a  constant 
near  the  origin.  If  T  is  some  measure  of  the  duration  of  e(t), 
then  as  T  increases,  a  point  is  reached  where  the  integral  in  the 
numerator  approaches  a  constant  value  (due  to  the  relative  narrowness 
of  Vnn(u)).  Meanwhile,  the  integral  of  the  denominator  grows  roughly 
as  T,  and  hence 


Lim  R0  -  Lim  -1 =  0  (5*+) 

T - >  CJ  I  —V  OO  1 


This  shows  that  as  the  duration  of  the  envelope  increases,  we  can 
expect  the  autocorrelation  function  of  s(t)  to  depart  percentage¬ 
wise  less  and  less  from  E  [/*ss('t'^  at  every  point. 

We  turn  next  to  the  spectrum  of  s(t),  where  it  is  found  that  such 
convenient  limiting  behavior  does  not  occur. 
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Chapter  IV 
NOISE  BURST  SPECTRA 


To  continue  this  investigation  of  the  properties  of  noise 
bursts,  it  is  now  appropriate  to  turn  to  a  consideration  of  the 
energy  spectrum  of  signals  of  this  type.  As  described  in  reference 
(a),  a  deterministic  transient  can  be  specified  by  its  Fourier 
integral: 


S(«)  =  — 2^' "  f  sCt)  e‘ 


(55) 


such  that 


s(t )  =  f°S(o)eiut  d» 


(56) 


Now  by  Parseval's  theorem  (reference  (a),  pp  33-39), 

2  n  y  s(«)  S CT  d»  =  J  s2(t)  dt 


(57) 


where  the  superscript  bar  represents  complex  conjugation.  Since  the 
right  hand  integral  is  simply  the  total  energy  of  s(t)  as  defined  in 
equation  (5),  the  expression 


Ct>ss(«)  =  2tt  S(u)  s(u)  -  2w  S(o)| 


(58) 


can  be  interpreted  as  an  energy  density  spectrum  indicating  how  the 
total  energy  of  s(t)  is  distributed  in  the  frequency  domain.  This 
is  to  say  that  if  s(t)  were  passed  through  an  ideal  low  pass  filter 
with  a  sharp  cut-off  frequency  uc,  then  the  total  energy  dissipated 
in  a  pure  one  ohm  resistance  following  the  filter  would  be 


E(«<  «c)  =  J 


d5 gs(«)du  s  2  /  <I>ss(«)dw 


(59) 


V  — 

Alternatively,  {I>cS(“)  can  be  interpreted  as  giving  the  total 
energy  dissipated  in  the  ubiquitous  one  ohm  resistor  by  a  band  of 
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frequencies  one  radian/sec  wide  centered  at  the  argument  radian 
frequency,  w  . 

For  a  signal  of  the  random  transient  type,  the  integral  of 
equation  (55)  can  be  written,  at  least  formally,  as  the  first  step 
of  obtaining  the  energy  density  spectrum  of  the  burst.  From  the 
same  considerations  that  led  to  the  formation  of  the  autocorrelation 
integral,  it  can  be  shown  that  the  integral  of  equation  (55)  exists 
for  all  members  of  the  ensemble  except  for  a  subset  of  probability 
zero.  But  since  every  member  of  the  ensemble  is  different  from  every 
other  member,  every  S(w)  will  be  different  also,  and  hence  there  is 
generated  an  ensemble  of  integrals  of  the  type  of  equation  (55)  de¬ 
fining  S(<*>)  as  a  new  random  variable.  Thus,  by  equation  (58), 

(t)ss(w)  will  also  be  a  random  variable  for  every  value  of  w  and  will 
not  possess  a  deterministic  value. 

The  next  step,  evidently,  is  to  compute  the  expected  value  of 
d>ss(M)*  By  definition. 


<j)ss(rj)  s  2tt  S(U)  s(u)  s  /  s(t)e*’^“tdt  /  s(x)e+^“x  dx 

/•  OO  p  o» 

S  -j*-  /  /  s(t)s(x)e“^“^t“x^dx  dt 


(60) 


-ga  —  o° 


Following  through, 


B 


(4>sg<»0  ■  S ir  I  f  e(t)e(x)  E  e'3“(t_x)dx  dt 

•  •  (61) 


Because  n(t)  has  been  assumed  ergodic  and  stationary, 


E  jn(t)n(x)J  s  ^(t-x)  =  ^(u) 


(62) 


where  u  *  t-x 
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Thus, 


©O  £>o 


[*s.(.U  =  -£r  f  [  e(t^e(t“u)  V'nnCu) 


e  ^wU  dudt 


/OO  z*00 

ywu>e"jw7e(t)e(t-u>dt  du 


(63) 


And  finally, 


••S&Ss<“D  =  -gt-  / 


e'3“u  du 


(6W 


which  can  be  interpreted  as  the  inverse  Fourier  transform  of 
E  pss(uQ  found  in  equation  (25).  Actually,  by  using  the  Wiener- 
Khinchin  relation,  it  is  possible  to  relate  the  expected  value  of 
<J>SS(<«0  to  the  energy  density  spectrum  of  e(t)  and  the  power  density 
spectrum  of  n(t).  Indeed,  with  ^n^^)  given  by  equation  (8),  the 
power  density  spectrum  of  n(t)  is  given  by  the  Fourier  transform 


*nr,(“>  = 


/o+ 

^nn^) 


e“J“r  dr 


(65) 


and  the  energy  density  spectrum  of  e(t)  is  given  by  the  transform* 


®ee<“>  =  27 


/  ?ee(r)e-3“r 


(66) 


The  inverse  transform  of  $ee(r)  is  given  by 


e-pp(r) 


■f 


<t>ee(»)eU*'  do 


(67) 


and  hence  equation  (61*)  can  be  written  as 


♦Similar  to  our  convention  in  Chapter  III,  the  upper  case  Greek 
letter^  (Psi)  will  be  used  to  denote  power  density  spectra,  whereas 
the  upper  case  letter  $  (Phi)  will  be  used  for  energy  density  spectra. 
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/C30 

<tee(“^3,’U  d.du 

-  -jU  (to-#) 

.  ~2^  /  j^nn^)  SW'Oe  J  d<r  du 


Now  setting  «'  *  <0  -  o-  ,  we  can  write  that 


E&ss<"3  =  Jv nn(u 


)  ^ee(w“  «0e“^“  u  dudw/ 


mO”  —  OO 


Using  equation  (65), 


-  /  <J>eek-w')d</ 


=  *rm(w)  »  ^ee^ 


(68) 


(69) 


(70) 


Thus,  the  expected  value  of  the  energy  spectrum  of  the  random 
transient  is  found  to  be  the  convolution  of  the  energy  spectrum  of 
the  envelope  waveform  and  the  power  spectrum  of  the  original  noise. 
At  least  as  far  as  expected  values  are  concerned,  it  has  been  shown 
that  one  can  here  apply  the  well  known  "folk  theorem"  that 
"multiplication  in  the  time  domain  corresponds  to  convolution  in  the 
frequency  domain." 

Continuing  with  this  Fourier  integral  approach,  it  is  possible 
to  investigate  the  variance  of  $ss(<*>)  at  every  point.  First,  by 
equation  (58), 


<S>  <.«.(«)  s  4tt^  S^(w)  S(oj) 


Now  using  equation  (55), 

»)pooo  om 

*ssM  S  i^5 JJJJ  s(t,s 


[(t)s(x)s(u)s(v)e"^w^t“x+J1“v^dtdxdndv 


(71) 


(72) 


o  •  o  o 


OOOOCWOO 


jjjj 


e(t)e(x)e(»i)e(v) 


o  o  0  o 


E[n(t)n(x)n(ju)n(vJe“^u^t"’x+p’"v^dtdxdpidv  (73) 
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Since  n(t)  is  assumed  to  be  Gaussian,  equation  (31)  still 
applies,  and  this  becomes 

O*  OO  (Xj  cxs 


B[**H  -  &Jfff  e(t)e(x)e(p)e(v) 


0000 


jVnn  ( t  -x )  Vnn  ( u-v )+  V^n  ( t-u)1^  x-v )+  %n(  t  -v )  %n(  x-uTj 


,-ju(t-x+^-v) 


dt  dx  dju  dv 


(7b) 


This  is  plainly  the  sum  of  three  integrals,  two  of  which  are  identical. 
The  first  and  third  are  of  the  form 


#»  #odo  oo 


II  =  13 


=  JL-  [.[[[$( t)e(x)e(p)e(v)V  (t-x)'^m(^-v)e“Ja(t“x+‘^“v) 

4"2  JJJJ  dtdxd*2dv 


0  C  o  o 


(75) 


Using  the  substitutions  a  -  t  -  x.  and  b  =  fi  -  v,  these  are  equiva¬ 
lent  to 

-1!  *  I3  * 


s  ^2  JVhn^a)e"*^ua  j  e(t  )e(t-a)dt  da  jv: nn05)e"^“b  |  e(Ai)e(;u-b) 

•-*  OO  O  *"  ^ 

s  £2  nWa^ee(a>e"3“a  4« J'f’m^eeibie'3^  ab 


cjndb 


y>(a)Vee(a)e~l»a  da 


,j  =|e  [a>ss(«)]| 


(76) 
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The  second  integral  is 


04*  OO  OO  00 


I2  *  ^ j'J'J'f'e(t)e(x)e(ii)e(v)  ^^(t-u)  ^nn(x-v)e 


0000 


00  00 


-j«(t-x+/u-v) 

dt  dx  d*i  dv 


(77) 


=  f  ^^^^(t^e^^^dt  du 

O  o 

&Q  CX* 

J  J e(x)e(v)Vnn(x-v)e+^“^x+v^dx  dv 

0  O 

The  first  integral  of  this  product  is  identifiable  as  the  complex 
conjugate  of  the  second,  so 

100  00  ,2 

J  J e(t)e(/u)  VWt“M)e~^t+^dt  dju 
Therefore,  from  equations  (7^),  (76),  and  (78), 


(78) 


B&>sS(“D  =  J1  +  x2  +  l3 


(79) 


*  2Je&5s(“D]2  +  1 5^  ^'e(t)e(n)Vnn(t-*i)e‘;!“(t+'j)at  ^ 


and  by  the  definition  of  the  variance, 

OO  OO 


Var[^s(^  =  jE&ss^^]|2  +  nn(t-h)e"^M^t+w^dt  dn 

0  0 


Var 


[4>ss<“>]  2  |e|*Ss(“)]) 


2 

(80) 


We  have  found,  then,  that  the  variance  of  the  energy  density 
spectrum  is,  for  each  «  ,  at  least  as  large  as  the  square  of  the  mean 
of  the  energy  density  spectrum  at  u  ,  or  alternatively,  that  the 
standard  deviation  of  <f>ss(&>)  is  always  larger  than  the  mean  -  regard¬ 
less  of  the  form  and  duration  of  e(t).  This  somewhat  surprising 
result  is  a  more  general  statement  of  the  relatively  well  known  fact 
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that  the  so-called  "periodogram"  which  can  be  drawn  from  an  empirical 
record  as  an  estimate  of  spectral  density  of  a  continuing  process 
does  not  converge  in  the  mean  to  the  true  spectral  density  as  the 
observation  period  increases  without  limit,  £see  reference  (b), 

pp  107-108). 

The  implication  of  these  findings  is  the  followings  If  a  large 
number  of  random  transient  sample  functions  are  taken  and  their 
energy  density  spectra  computed  or  measured,  the  arithmetic  average 
of  <J>gg(w)  will  tend  toward  the  expressions  of  equations  (64)  or  (70) 
for  every  to  .  Because  of  the  magnitude  of  the  variance,  however,  it 
is  quite  likely  that  the  value  of  $ss(“)  for  a  given  sample  will  be 
nowhere  near  the  mean.  Consequently,  at  a  given  value  of  frequency, 
there  will  in  general  be  a  wide  variation  of  the  energy  density  spec¬ 
trum  from  sample  to  sample.  Unlike  the  autocorrelation  function 
which  becomes  percentagewise  more  precise  as  the  duration  of  e(t) 
increases,  the  empirical  spectrum  of  each  sample  function  remains 
widely  distributed  regardless  of  the  form  of  the  envelope. 

In  Appendix  B,  an  alternative  approach  to  the  noise  burst 
spectrum  is  based  on  the  assumption  that  the  Wiener -Khint chine  rela¬ 
tion  applies  to  signals  of  this  type.  The  results  obtained  are 
identical  to  those  presented  above,  and  it  appears  that  such  an 
assumption  Is  valid. 
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Chapter  V 

EXAMPLES  OF  NOISE  BURST.  AUTOCORRELATIONS  .AND  SPECTRA 


To  demonstrate  the  application  of  the  principles  derived  above, 
this  section  will,  be  devoted  to  the  presentation  of  several  examples 
of  the  calculation  of  noise  burst  autocorrelation  functions  and  energy 
spectra.  For  this  purpose,  two  idealized  noise  spectra  will  be  treated: 
fhe  rectangular  low-pass  broad  band  spectrum,  and  the  rectangular 
narrow  band  spectrum.  These  will  each  be  modulated  by  two  envelope 
functions:  a  rectangular  pulse  of  duration  T  seconds;  and  an  expo¬ 
nential  decay  with  a  time  constant  of  T  seconds.  This  yields  four 
cases: 

A,  Broad  Band  Spectrum  with  Rectangular  Envelope 

B.  Broad  Band  Spectrum  with  Exponential  Envelope 

C,  Narrow  Band  Spectrum  with  Rectangular  Envelope 

D.  Narrow  Band  Spectrum  with  Exponential  Envelope 


The  two  noise  spectra  and  their  associated  autocorrelation  funcr 
tions  are  shown  in  Figures  2  and  3,  For  the  rectangular  low  pass 
broad  band  spectrum, 


*nnM  =  No,  -«2^w2 
=  0,  elsewhere 
By  the  Wiener-Khintchine  relation, 


Thus, 


sin  woT 

*nn(T)  2N0  “2  "  TTTr ~~ 


(81) 


(82) 


(83) 


For  the  rectangular  narrow  band  spectrum, 

=  No»  wl-“  -  “2 


~  No*  -w2  s  oo  s  ~  ojj_ 
=  0,  elsewhere 


(84) 
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Thus 


fnni 


r  102 

r)=2/  K 

Jus-, 


o  cos  cor  d  w 


(85) 


=  4  A  a)  No 


sin  a  cor 
Ac or 


t 


where 


A  CO  = 


Ob  = 


<02-  wi 
2 

^2  "Hoi 


(interpreted  as  the  half 
bandwidth) 


(interpreted  as  the  center 
frequency) 


(86a) 


(86b) 


Turning  now  to  the  first  of  the  envelope  functions,  consider  the 
rectangular  pulse  shown  in  Figure  4. 


e(t)  =1,  OstsT 
=  0,  elsewhere 

The  autocorrelation  function  is  computed  from  equation  (4) 


(87) 


<Pee(r)  =  J e(t)e(t  +t)  dt 


=  T  -  |t|  » 


-T  *  r  *  T 


The  energy  spectral  density  can  be  found  as  either  the  Fourier 
transform  of  equation  (88)  or  as  2ir  |E(co)|2  where  E(co)  is  the 
Fourier  transform  of  the  envelope  waveform. 


<x>eew  =  -il  iiaLfl'k 
2  T  ( “V2)2 


(89) 


One  more  quantity  is  needed  in  the  study  of  rectangular  noise 
bursts,  and  this  is  the  function  V’pptu, r)  defined  by  equations  (40) 
and  (41).  If  r  is  positive, 

p(t,r)  =  e(t)e(t  +r)  =  1,  0<t<T-7  (90 


Therefore, 


Vu’ r  > 


■/ 


T-r-  u 


p(t,  T  )p(t-u,  T  )dt 


/I  « 

1  dt,  pt  0 
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^pp  (u,  T)  =  T  -T-  u  fox  rand  u>0 

O^uST  -  r 

Since  (u,  T)  is  an  even  function  of  both  rand  u, 

^pp(u,  r)  =  T  -  Ir|  -  |u|  ,  [|r|+|u|]<T 

=  0,  elsewhere 


(91) 


(92) 


Let  the  exponential  envelope  function  portrayed  in  Figure  5  be 
defined  as 

-t/_ 

e(t)  =  e  0  ^  t  <  00  (93) 

=  0,  elsewhere 


Defining  the  autocorrelation  function  precisely  as  in  the  case 
of  the  rectangular  pulse,  we  find 


%AV) 


,  -oo<V<< 


(94) 


The  energy  density  spectrum  becomes 


Oee(co)  = 


_  T2 


27f  co2  T2  +1 


(95) 


and  finally. 


^ppfu*  T1 )  = 


_  T 


-T-(|ui+  in  ) 


,  -oo<u,r 


(96) 


For  convenience,  these  functions  are  collected  in  Table  I* 
AUTOCORRELATION  CALCULATIONS 

1.  RECTANGULAR  PULSE  OF  BROAD  BAND  NOISE .  By  equation  (25), 

-  %nW  %e(T) 

and  for  this  case,  the  autocorrelation  functions  for  the  rectangular 
pulse  and  the  broad  band  noise  spectrum  are  found  in  Table  I.  Using 
these  expressions, 

E  [%s(r  )]  =  2N0  to2  (T  -jri)  ,  -T^r  <T  (97) 

=  0,  elsewhere 
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FIGURE  5. 


EXPONENTIAL  ENVELOPE  AND  ASSOCIATED  AUTOCORRELATION 
FUNCTIONS. 
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From  equation  (43)  can  be  found  an  expression  for  the  variance  of 

»SSW: 

VarjU  (-r)j  =  2  f  f  pp(u,r)  [y.nn2(u)+^nn(u+r)  fnn(u-r)]  du 
J  o 

Referring  again  to  Table  I, 

(98) 

r  r T-  H  r  2  1 

VarUss(T)l  =  8N02«22  /  (T-T-u)  sin  “2  »  .  sin^lr+u)  sin^jr-u) 

L  J  Jq  w22  p2  w2 (r+u )  w2lr-u)“ 

The  standard  deviation  of  v’ssfr)  is,  of  course,  the  positive  square 
root  of  this  expression. 

For  more  generality  of  presentation,  the  expressions  of  equa¬ 
tions  (97)  and  (98)  will  now  be  normalized  and  parametized.  We  note 
first  that  since  ¥>Ss(Oj  is  equal  to  the  total  energy  of  the  noise 
burst  (as  defined  in  equation  (27)),  eTv’ss  (0)1  can  be  described  as  the 
average  total  energy:  L 

ETav  =  E  Us  (0)]  (99) 


For  the  present  case, 


^Tav  =  2No  <*>2T 


(100) 


which  makes  sense  since  it  turns  out  to  be  equal  to  the  average  power 
per  radian  times  the  bandwidth  times  the  duration  of  the  transient. 
ECv’ss(r)^  will  be  normalized  by  dividing  by  Ejav  and  Var [Uss (t)3 
by  dividing  by  (Ejav)2»  This  latter  step  puts  the  standard  deviation 
of  v»ss(r)  in  the  same  units  as  the  mean.  Thus: 


En  [  3  (T  j]  =  (1  -  — — )  -T  <  T  <  T 

w2  T 


(101) 


,  elsewhere 


t-  in 

nUssCDl  =  2/  <+-  h  -  J-> 

L  J  J  11  i*  [_w2.  U2  "2(T+  p)  "2.1  T-zu) 


(102) 
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Now  let  the  following  parameters  be  defined: 


(103a) 


p  -  -21-  t  i,ev  the  argument  as  a  fraction  of  the  total  pulse 
T  length 


q  = 


u>2 


o>2T 

_ ,  or  the  number  of  periods  of  the  highest 

2  *  noise  frequency  contained  in  a  pulse 

length.  (103b) 


The  normalized,  parametized  equations  can  now  be  written  as 


Epj^ssIP}] 


(1-p) 


sin  2  7r  pq 
2tt  pq 


(104) 


Var, 


rl-P 

[*ss(p)]  =  2  J(l-p-x) 


sin^  2?rqx 
(27rqxj^ 


sin  2irq(x-fp)  sin  2Trq(x-p) 
2irq  (x+p)  2irq(x-p) 


(105) 

dx 


By  equation  (52),  the  ratio  of  the  standard  deviation  atT=  0  to 
the  expected  value  of  ^ss(t)  at  7=  0  is 


(106) 

and  it  should  be  recalled  that  this  ratio  is  the  maximum  value  of  the 
standard  deviation  expressed  as  a  fraction  of  the  iriean.  It  should 
be  apparent  from  the  foregoing  normalization  and  parameterization  that 


-x)  sin22.J-qx-.  dx  (107) 

(27rqx  )2 

The  functions  of  equations  (104),  (105),  and  (107)  have  been 
programmed  for  evaluation  on  the  IBM  7090  computer  at  NOL,  the. 
latter  two  expressions  requiring  the  use  of  numerical  integration 
subroutines.  Computer  generated  plots  of  the  normalized  expected 
value  of  vss(p)  are  found  in  Figures  6  and  7  for  this  case  when 
q  =  0.1,  1.0,  10.0,  and  100.0.  It  should  be  realized  that  all  such 
autocorrelation  functions  are  even  functions  of  p  and  disappear 
for  |  p  |  >1.  Also  note  the  change  of  scale  for  the  p  axis  when  q=100 
(in  Figure  7).  Figure  8  shows  a  typical  plot  obtained  by  graphing 
on  the  same  axes  Bp  fyss(p)]  and  this  mean  value  plus  and  minus  the 
standard  deviation  of  v>ss(p")  as  calculated  as  the  square  root  of 
equation  (105).  In  some  sense,  such  a  graph  provides  a  rough  idea 
of  how  closely  the  distribution  of^>Ss(p)  cleaves  to  the  mean  as  a 


R0  =  Ro(q) 
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NORMALIZED  AUTOCORRELATION 
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RECTANGULAR  BROAD  BAND  Q  =  10 


FIGURE  8.  RECTANGULAR  PULSE  OF  BROAD  BAND  NOISE 

AVERAGE  AUTOCORRELATION  t  AUTOCORRELATION 
STANDARD  DEVIATION 
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function  of  p.  Generally,  the  standard  deviation  will  be  largest  for 
p  s  o  and  then  decrease  as  p  increases.  This  decrease,  however,  is 
not  always  monotonic,  and  the  expression  R0(q)  remains  perhaps  the 
best  single  estimate  of  the  tightness  of  the  distribution  of  <fSs(p)« 

A  plot  of  R0(q)  is  portrayed  in  Figure  9  and  mirrors  the  result 
predicted  in  equation  (54).  Since  the  parameter  q  is  a  measure  of 
the  length  of  the  envelope  pulse  in  terms  of  the  period  of  the  highest 
noise  frequency,  it  will  be  directly  proportional  to  the  envelope 
"duration"  discussed  in  the  derivations  of  equation  (54).  Thus  as  q 
increases,  we  can  expect  R0  to  decrease  sharply,  and  indeed  this  is 
what  has  been  found  here.  As  the  pulse  length  increases  relative  to 
the  period  of  the  highest  noise  frequency  present,  the  distribution 
of  ^fss(p)  clusters  closer  and  closer  to  the  mean  for  all  p,  and  the 
value  of  the  mean  at  p  becomes  a  better  and  better  estimate  of  the 
actual  value  of  <fss(p). 


2.  EXPONENTIAL  PULSE  OF  BROAD  BAND  NOISE.  Again  using  equa¬ 
tions  (25)  and  (43)  and  the  appropriate" 'entries  from  Table  I,  it  is 

hat 

=  N0^Te  S;in<^r  (108) 

w2r 


found  immediately  t 
E 

and  that 


Var  [Vss(r)]  =  2  N02^22  t  e~2^T 


(109) 


oo 

I/' 


-2u/t 


sin2^  sin  ^2(u+r)  sin  u>2(u_r)' 
L^2  »4  +  w,(u+r) 


It  follows  directly  that 


ETav  =  E  [<».«»]  =  N0 T  =  N0^  fe 2(t)  dt 

Jo 

As  before,  the  functions  of  equations  (108)  and  (109)  will  be 
normalized  and  parametized.  Defining, 


du 


(110) 


and 


rt> 

p  a  — =~  ,  i.e.,  the  argument  as  a  fraction  of  the  envelope 

time  constant  (111a) 

q  =  — —  ,  the  number  of  cycles  of  frequency  ^  contained 
2  7T  in  the  time  constant  T  (111b) 


the  functions  become 


=  e-p  sin  2<*pq 
2?f  pq 


(112) 
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RECTANGULAR  BROAD  BAND 
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FIGURE  9.  RECTANGULAR  PULSE  OF  BROAD  BAND  NOISE 

NORMALIZED  AUTOCORRELATION  STANDARD  DEVIATION. 
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r°° 

fpss(p)]  = 2e_2p/ e_2x  |j 


(113) 


sin^  27rqx  +  sin  2*q(x- 
{2Wqx)2  27rq(x+p) 


in  2*g(x-p)  , 
2  ft q(x-p) 

mi 


and  following  through, 


R0tq)  =  2 


1  sin^^qx 
(2Tqx )2 


(114) 


Equations  (112)  and  (114)  have  also  been  programmed  on  the 
IBM  7090  computer,  and  the  results  appear  in  Figures  10,  11,  and  12. 
As  before,  the  average  autocorrelations  are  presented  for  q  =  0.1, 
1.0,  10.,  and  100.  Note  again  that  as  q  increases,  R0  falls  toward 
zero.  As  the  time  constant  of  the  envelope  increases,  the  distribu¬ 
tion  of  9? ss  (p )  closes  more  tightly  about  the-  mean,  as  one  would 
expect  from  previous  considerations. 

3.  RECTANGULAR  PULSE  OF  NARROW  BAND  NOISE.  Following  exactly 
the  same  procedure  as  before  and  again  referring  to  Table  I,  one  can 
write  that 


e  «ess(r)  a  4  No^T-lTl  ) 


r  i  fT~lrl  r  . 

Varj?ss(r)J  *  32  N o\u2 J  (T  -lt|-u)  cos2o>oa 


a^t 


cos  o)0  T 


(115) 


(116) 


+  sinAa)(M-T)  sinA^+r) 

(m-7)  Aw(^+r) 


where  as  before 


^o  = 


Now  from  equation  (115), 


cos^(4-r)cosd0(M+r)  d<u 


AcJ  = 


cu2 


ETav  =  4  N0  &ojT 

which  is  what  one  would  expect  intuitively.  For  this  case,  the 
following  parameters  will  be  defined: 


(117) 
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EXPONENTIAL  BROAD  BAND  Q  =  0.1000 


P  (TAUAIME  CONSTANT) 


EXPONENTIAL  BROAD  BAND  Q  =  I 


P  (TAUAIME  CONSTANT) 


FIGURE  10.  EXPONENTIAL  PULSE  OF  BROAD  BAND  NOISE 
AVERAGE  AUTOCORRELATION 
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EXPONENTIAL  BROAD  BAND  Q  =  10 
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P  (TAUAIME  CONSTANT) 


EXPONENTIAL  BROAD  BAND  Q  =  100 
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FIGURE  11.  EXPONENTIAL  PULSE  OF  BROAD  BAND  NOISE 
AVERAGE  AUTOCORRELATION 
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EXPONENTIAL  BROAD  BAND 


2.0 


1.5 


1. 


0. 


0.1  1  10  100  1000 
Q  (CYCLESAlME  CONSTANT) 


FIGURE  12.  EXPONENTIAL  PULSE  OF  BROAD  BAND  NOISE 

NORMALIZED  AUTOCORRELATION  STANDARD  DEVIATION 
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P  -  T  »  the  argument  as  a  fraction  of  the  pulse  length  (118a) 

q  =  7T  =  *  Jhe  nu7nber  of  cycles  of  the  center 

2  r  frequency  round  in  a  pulse  length 

Z  =  izf*  the  ratio  of  bandwidth  to  center  frequency 


(118b) 

(118c) 


^rrofCfnd^iae^at510"5’  “  “  f°Und  for  a  ^^tangular  burst  of 


Ep^Jp)]  =  (i-p)  si!H35B  tos2,pq 


(119) 


r  I 

r  X-P  r 

jq’ssfp)  J= 

2  (1-P-x) 

•'o  ^ 

sin  tf  qzx  o 

cos  27rq* 


ro(9iZ)  =  2 


4.  Iinrqz(x-p)sin7rqz(x+p)  _  ,  1 

^rqzU-p)  ~T^z{x+p)  C°S  2rq^x"P)cos  2vq(x+p) I  dx 


sin27rqzx  _ .  2  _  . 

W^C0S  2rqxdx 


Typical  results  when  equation  (119) is  pvaina+oH  , 

computer  appear  in  Fiaures  anrf  1/  r  1S  evaluated  on  a  digital 

Ro  <q,z)  is  plottad  af a  fu^tiCn  of'o  vf+b  fh  Z  =  °-1*  1,1  Fi9ure  1EV 
z-value  of  2?3,  which  may  seem  like  an  £d  choice  The 

So?ee^a°?dR  (o  z)  aCtaVe  ^and  sPe^rum,1a  S^Jl/bP Jelified’ 

£*%&%&  (dq4c2ii^P?rS^e5lyr°m^qraiSdre|Ir tL°VU  *  2 

z,  corresponding  to  wider  and  wider  band?widths.th  1  961  ValU6S  °f 

usin9  tha  propar 


-  r 

[<?ss(r)]  =  T 


2N0T  Aa>e  1  sin  ACJ-r 


wTiere  A  cu  and  cd0  are  defined  as  before 


aw  r 


cos  oL  T 
0  1 


(122) 
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RECTANGULAR  NARROW  BAND  Q  =  0.1000 
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RECTANGULAR  NARROW  BAND  Q  =  1 


FIGURE  13.  RECTANGULAR  PULSE  OF  NARROW  BAND  NOISE 
AVERAGE  AUTOCORRELATION  FOR  Z  =  0.1 
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RECTANGULAR  NARROW  BAND  Q  =  10 
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Var[<pss(r)]=  8N02Tao;: 


~2r/jJe  "2u/T 


s,in2AO)u.  cos2 
(Aa)u)2 


(123) 


sinAatfu-r)  sin Atofu+r) 


Au)  ( u-  r) 


AU)  (U+77 


cosotiu-r)  cos  at  (u+r) 


Ejav  =  2NqT  A  cO 


(124) 


The  parameters  for  this  case  are  precisely  the  same  as  those 
defined  for  the  rectangular  narrow  band  ’'urst  with  the  exception 
that  T  is  now  to  be  taken  as  the  envelop  time  constant*  This  yields 


)[«?ss(p)]  = 


7rqzp 


cos  2?rpq 


[•..H  ■  -  te 


—  cos2  2irqx 

)2 


(125) 


(126) 


4  -sinyqz fozP-)  sinyqzfx+pj.  COs  27rq(x-p)cos  2flq(x+p)  dx 
rqz(x-p)  Tt  qz(x+p) 


Ro(q»z)  =  -[I  Is  cos2  2^rqx  dx 

U  (TTqzx  )2 


(127) 


Typical  plots  of  these  equations  are  found  in  Figures  16f  17, 
and  18*  They  tend  to  bear  out  the  generalizations  previously  noted 
for  the  other  cases. 


SPECTRUM  CALCULATIONS 

The  fexpected  value  of  the  energy  density  spectrum  of  a  noise 
burst  for  some  argument  a;  has  been  found  to  be 

E  [*„(*)]  =  f*nn((o)  dto' 

—CD 

=  <8>  'M"5 

by  equation  (70)  and  is  interpreted  as  the  convolution  of  the  power 
spectrum  of  the  noise  and  the  energy  spectrum  of  the  envelope 

b-7 
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NORMALIZED  AUTOCORRELATION  NORMALIZED 


FIGURE  17  .  EXPONENTIAL  PULSE  OF  NARROW  BAND  NOISE 
AVERAGE  AUTOCORRELATION  FOR  Z  =  0.1 
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Q  (CYCLESAIME  CONSTANT) 


FIGURE  18.  EXPONENTIAL  PULSE  OF  NARROW  BAND  NOISE 

NORMALIZED  AUTOCORRELATION  STANDARD  DEVIATION 
PARAMETERIZED  BY  BANDWIDTH. 
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waveform.  Since  both  members  of  the  convolution  are  even  functions, 
the  expected  value  of  3>ss(w)  can  be  written  as  the  cross-correlation 
of  *nnM  and  $ee(<d): 


f 60 

[(tssM]  =  *  nn 

~cu 


(eo*  )  c i>ee[(j0'-u))  du;1 


(128) 


Obtaining  the  spectral  functions  of  the  various  envelopes  arid  noises 
from  Table  I.  it  is  straightforward,  although  often  tedious  to 
calculate  E  [<t>ss(^)]  for  the  examples  treated  here. 

1.  RECTANGULAR  BURST  OF  BROAD  BAND  NOISE.  Using  formula  (128), 


K«1  •  ¥[ 


U3 


NoT2  j  + 

(^2 


sin2  co*  V2 

(w*  V2)2 


dw» 


(129) 


Now  introducing  the  same  parameters  defined  in  the  study  of  the 
autocorrelation  function  for  this  case  (equations  (103))  and  adding 
another: 


r  =  -Eg-,  the  argument  as  a  fraction  of  the  highest  noise 
frequency,  a  normalized  and  parametized  version  of  equation  (130) 
can  be  written  as 


1 


f  q(r+l)  . 

1  I  sin-i 

2  -'q(r-l)  M)2 


(131) 


such  that 


(132) 


From  the  normalized  expression,  one  obtains  the  spectral  density 
for  a  one  radian  band  aico=  r«2  by  writing 


(133) 


=  N0T 


r  q(x+l) 
q(x-l) 


s-4agB  dx 

(r  x)2 
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Now  integrating  equation  (131),  it  emerges  that 

E[<  <r>]  = 


(134) 


1  )  1-r  sin  2?rqr  sin  2rtq-2  cos  2-aqr  cos  2ttq 

2^2  j  5(3£l} 


+  ^f  {Si  [2'7>'c>tr+1)]  -Si[2rq(r-1)]  } 
where  Si  is  the  sine  integral  function  defined  as 


Si(x)s  /X 
Jo 


dx 

x 


(135) 


(See  reference  (d)) 


Evidently,  if  the  spectral  width  of  the  envelope  is  narrow  com¬ 
pared  to  that  of  the  original  noise,  the  convolution  of  equation  (128) 
will  just  return  an  approximation  to  the  latter.  Such  a  condition 
implies  q  verv  large  (hy  reciprocal  spreading  arguments),  and  in 
equation  (134),  the  first  term  will  become  negligible  in  comparison 
with  the  second.  Since  for  large  positive  argument,  Si(x)  «  7r/0, 
and  for  large  negative  argument,  Si(x)«  -  * 


[«ks<r>l  -2V  jSi  [2^(r+1)]  -Si  £ 2?rq(r-l) j  | 


(136) 


«  1/2,  -1  <r  <+l 

0,  elsewhere 

which  checks  with  the  intuitive  solution. 


2.  EXPONENTIAL  BURST  OF  BROAD  BAND  NOISE.  Using  formula  (128), 


N0T *  [“>?■+“>  dw- 


E  L*ss M  =  -j 


Uf  ^  +1 


(137) 


where  T  is  the  decay  time  constant.  Normalizing  and  parasitizing 
as  before, 


[*>>]  ■  / 


q(r+l) 

dx 

4T2  x2+l 


(138) 


q(r-i) 


where  r  = 


?2 
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Integrating, 


[*»<r)] = 


.1  + 


47Tq 

47f2  q2  (j2 


-1). 


(139) 


3.  RECTANGULAR  BURST  OF  NARROW  BAND  NOISE,  The  condition 
integral  can  be  written  as 


-D0+a>+  &a) 


5ln2£  Tf-  **  * 

L  J  27r  J  («•  T/2  j2 

-ft)0  +0)- 


o  ru)+^)  +  ^UJ 

NoT  f  sin^co'  T/*> 

2^  J  U»  T/2)2 


(140) 


W0+W-,ACU 


Defining  as  before, 

_  <a0T 

q  -  2?r  » 

this  can  be  written  as 


2 


[*«(r)]  ‘TzJ 


q(r-l+  z/2) 
f  sin2  ttx 

■J  (TTx)2 

q(r-l  -z/2) 


+  2'/, 


q(r+l+  z/2) 

sipi^.x  dx 


/  (^x) 

q  ( r+1—  z/2) 


(141) 


in  the  parametized  unit  energy  version.  In  this  fox-m, 

/GO 

E  [*.,(*)]  dr  =  1,  and  the  relationship  between  E J<5>ss  (c*>)j  and 

—  CO 

E^ss(r)]  can  be  expressed  as 

=  2  N0Tz  E  [<t>ss(*  =  ^  )] 

=  £lsv  H  ks(r^A  )j  (142) 

UJn  L  6>o  J 


At  any  r-ate,  integration  yields 
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Fcj)  (r)l  ss  1 

1-* 

t 

1 

1 

L  SS  J  8  7f2 

1 

;(x-i)2  -  z2^; 

q 

(x+l)2  -  JE2/  1 

L  '4 1 

1  f,2(x-l)sin  2lfrq(x-l)  sin  'Wqz  +  z  cos  2  nr q(r-l)cos  rtqz 
4ir2qz  \  [(r-1  f  -  ,8/4  ]  ' 

+  2(x+l)sin  -2y  q(x+l)  sin  ^qz  +  z  cos  27t,q(r-fl)  cos  rfqz  I 

[(x+1)2  -  z2/  j  / 


|si  [2wq(l+x+z^)]  -  Si  [2 irq  (1+x  -  z/2  )jl 


It  can  be  shown  that  fox'  q  laxge, 

E'[*”(r,r&'  •’  -1  ■  z/2£xs-^  +  V2 


(143) 


(144) 


IT'  i-vas*si-*«/s 

i  0,  elsewhexe 


as  it  should. 


4.  EXPONENTIAL  BURST  OF  NARROW.  BAND  NOISE.  Evidently, 


-CO  4.  CO  4.  A  to 
0  T  T  “ 


E[‘ 


H  4>ssM  ■« 


1  _  No*2  J  M  N0T2 

J  ’  2*  J~u J2  T2  +1  +  ”2^ 


+  <o>_ 


too+cO+ a<*> 


du/ 

"72-2" 
o>  T  +1 


%  t^-Aco 


(145) 


54 


-  o 


:#•  *  * 


NOLTR  67-25 


Proceeding  as  before, 


/’q(r-1+2/2)  /.q(r+l«/2 ) 

istnj  -i  /-  —  +  i  / _ 

J  z/^2x 2*i  2  /  4^77 

•'q(r+l-y2  ) 


(H6) 


such  that 


[®ssM]  =  E  <i-)J=  No  Tz  E  [<(,)] 


(147) 


Integration  of  equation  (146)  yields 


=  — 
L  J  2  : 


2/?l'qz 


1  +  4  ir2  q2£(r-l)2  -z2^  j 


+  tan' 


27fqz 


1  +  4?r2  q2[(r+l)2  -  z2/4] 


(148) 


ms.**?*.. 


This  concludes  our  consideration  of  illustrative  examples. 
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Chapter  VI 

CONCLUSIONo  RBCOltfESpATipHS 


The  foregoing  discussions  have  attempted  ic  „f  •v»*mulat.e  a  means 
for  obtaining  statistical  descriptions  of  the  proporv-oc  of  a  class 
of  random  transients  perhaps  well  described  as  "noise,  bursts"*  In 
so  doing,  it  has  been  necessary  to  thread  a  careful  course  between 
the  methodology  that  applies  to  deterministic  transients  on  the  one 
hand  and  that  intended  for  continuing  stochastic  processes  on  the 
other.  Attention  has  been  restricted  to  the  time  average  auto¬ 
correlation  function  and,  the  energy  density  spectrum  for  given  values 
of  time  displacement  and  radian  frequency,  respectively. 

Given  a  random  transient  drawn  from  the  ensemble  of  all  such 
signals  available  from  the  generating  mechanism,  it  is  possible,  at 
least  formally,  to  compute  the  time  average  autocorrelation  function 
by  the  familiar. process  of  displacement,  multiplication,  ;and  integra¬ 
tion.  It  is  similarly  possible  to  calculate  the  energy  spectral 
density  of  the  sample  function  by  either  Fourier  transformation  of 
the  measured  autocorrelation  function  or  by  a  Fourier  integral  treat¬ 
ment  of  the  function  itself.  Now  since^  each  of  the  transient  sample 
functions  is  different,  it  is  hardly  surprising  to  find  that  each  of 
the  measured  autocorrelations  and  spectra  will  be  different  also. 

This  implies  that  these  latter  functions  are,  for  every  value  of 
their  arguments,  random  variables,  in  the  sense  that  we  lack  exact 
a  priori  knowledge  of  their  values,  and  therefore  cannot  predict  the 
autocorrelation  function  and  spectrum  of  each  individual  transient 
with  exactitude.  Thus,  the  autocorrelation  function  and  spectrum  of 
a  fandom  transient  must  be  described  by  probability  distributions 
parametized,  in'  a  sense,  by  the  arguments  of  the  functions.  This 
investigation  has  not  attempted  the  derivation  of  the  form  of  these 
distributions  at  each  point,  but  has  been  restricted  to  a  calculation 
of  the  means  and  variances  of  the  spectra  and  autocorrelations  as 
functions  of  the  arguments. 

If  the  random  transients  treated  here  are  modeled  as  the  product 
of  an  envelope  waveform  and  a  continuing  random  process,  the  calcula¬ 
tion  of  the  means  and  variances  described  above  is  straightforward. 

In  the  case  of  the  resultant  autocorrelation  function,  the  mean  value 
at  every  point  is  found  to  be  the  product  of  the  autocorrelation 
functions  for  the  envelope  and  the  original  noise  process.  This 
result  agrees  with  intuition  and  is  similar  to  that  found  in  seeking 
the  autocorrelation  function  of  the  product  of  two  independent  random 
processes.  The  variance  of  the  resultant  autocorrelation  function 
can  similarly  be  expressed  in  terms  of  the  autocorrelation  functions 


56 


NOLTR  67-25 


of  the  original  signals,  but  admits  no  ready  intuitive  explanation. 

The  most  important  characteristics  of  the  variance,  however,  ure 
first  that  it  achieves  its  maximum  value  at  the  origin  (where  the 
mean  corresponds  to  the  total  average  energy),  and  second  that  as 
the  suitably  defined  "duration"  of  the  transient  increases  with 
respect  to  a  typical  noise  period,  the  standard  deviation  becomes  a 
smaller  and  smaller  percentage  of  the  mean.  This  merely  reflects  the 
fact  that  as  the  transient  lengthens,  the  signal  looks  more  and  more 
like  a  continuing  random  process,  for  which  the  autocorrelation 
function  <?Y  every  point  is  well  defined. 

Turning  to  the  mean  value  of  the  spectrum  at  a  point,  a  satisfy¬ 
ing  and  intuitive  result  is  found.  Since  the  random  transient  is 
formed  by  the  multiplication  of  two1  signals  in  the  time  domain,  one 
might  expect  the  resulting  energy  density  spectrum  to  emerge  as  the 
convolution  of  the  two  corresponding  spectra  in  the  frequency  domain. 
When  speaking  of  the  mean  value  at  each  point,  this  is  found  to  be 
the  case.  The  variance  of  the  spectrum  has  also  been  treated  and 
expressed  in  terms  of  the  parameters  of  the  original  signals.  In 
contrast  to  the  autocorrelation  standard  deviation,  that  of  the  energy 
density  spectrum  is  always  larger  than  the  mean  at  a  point,  regardless 
of  the  duration  of  the  transient.  Thus,,  the  measured  spectrum  does 
not  converge  in  the  mean  to  the  value  predicted  by  the  convolution* as 
the  transient  lengthens,  and  it  appears  that,  at  least  on  the  basis 
of  a  pointwise  comparison,  a  large  spread  of  measured  spectra  will 
always  'be  observed.  This  result,  as  was  pointed  out  previously,  is 
the  major  ‘defect  of  so-called  periodograra  analysis,  which'  can  be 
treated  as  a  special  case  of  the.  problem  faced  here. 

It  has  been  stressed  throughout  this  report  that  the  results 
derived  apply  only  at  specified  points  on  the  autocorrelation  and 
spectral  functions  wftep  no  knowledge  is  assumed  about  the  behavior  of 
the  function  at  other  points.  In  other  words,  the  means  and  variances 
derived  here  stem  from  unconditional  probability  distributions  for 
every  argument  value,  in  which  the  behavior  at  a  point  is  treated  in 
isolation.  For  this  reason,  it  is  risky  to  attempt  to  extend  the 
present  findings  to  describe  the  extent  to  which  the  empirical  auto¬ 
correlations  and  spectra  §£  a  whole  are  predicted  by  the  calculated 
mean  values.  One  could  envisage,  for  example,  a  sample  function  that 
yielded  an  empirical  autocorrelation  function  quite  similar  in  form 
to  the  expected  value  but  having  one  or  two  pathological  points  of 
substantial  disagreement.  The  examination  of  this  sort  o*f  effect 
requires  the  study  of  conditional  distributions  of  the  autocorrelation 
functions  and  spectra,  or  alternatively,  the  determination  of  the 
joint  density  functions  of  their  values  at  two  or  more  arguments. 

This  has  not  been  done  here  and  remains  a  large  and  interesting  area 
for  future  , invest igation*  At  present,  we  must  limit  ourselves  to  the 
consideration  of  single  points  and  resist  the  temptation  to  extend 
the  pointwise  conclusions  to  the  autocorrelations  and  spectra  in  their 
entirety". 

The  problem  of  power  spectral  estimation  from  empirical  records  is 
an  area  that  resembles,  in  many  ways,  the  study  of  random  transients. 
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In  both  fields  one  must  work  with  finite  length  segments  of  random 
processes  whose  good,  behavior  arises  primarily  from  their  extending 
in  time  from  -  ©o  to  +oo  and  this  leads  to  computational  difficult¬ 
ies,  A  good  many  of  the  techniques  of  power  spectrum  measurement  can 
probably  be  applied  to  the  present  st  dy.  Blackman  and  Tukey  (refer¬ 
ence  (f)),  for  instance,  study  the  problem  of  joint  estimation 'of 
neighboring  points  on  an  empirical  power  spectrum  and  derive  express¬ 
ions  for  smoothed  spectral  density  estimates  which  abandon  the  concept 
of  point  estimation  in  favor  of  band-wise  calculations  which  display 
a  higher  statistical  reliability.  This  appears  to  be  a  particularly 
fruitful  approach  for  the  class  of  problems  treated  here  and  may  lead 
to  more  meaningful  prediction  of  the  spectrum  and  autocorrelation  of 
a  sample  noise  burst. 

A  closely  related  area  is  the  derivation  of  a  linear  system  theory 
for  signals  of  this  type.  In  effect,  this  would  indicate  the  results 
to  be  expected  when  random  transients  are  subjected  to  filtering  and 
other  linear  operations.  It  would  relate  directly  to  practical 
problems  of  measurement,  detection,  and  interference  elimination,  and 
may  even  lead  to  the  development  of  optimum  linear  filters  in  the 
Wiener  sense.  These,  are  only  a  few  of  the  new  directions  that  can  be 
followed  in  further  work  on  noise-burst -like  waveforms. 

Despite  the  limitations  set  forth  above,  the  present  theory  has 
several  interesting  implications  for  the  processing  -and  measurement 
of  random  transient  signals.  It  indicates  to  some  extent,  for  example, 
thie  degree  of  reliability  that  can  be  assumed  in  basing  a  measurement 
program  upon  a  given  number  of  sample  functions.  The  average  auto¬ 
correlation  function  for  a  given  displacement  or  the  average  spectral 
density  for  a  given  value  of  w  can  be  found  by  averaging  a  sufficiently 
large  number  of  empirical  calculations.  The  variance  of  these  averages 
can  in  turn  be  estimated  by  turning  to  the  theory  set  out  here. 
Evidently,  the  longer  the  transient,  the  fewer  the  sample  functions 
required  to  give  reasonably  good  knowledge  of  the  average  autocorrela¬ 
tion  at  a  point.  This  consideration  does  not,  however,  apply  to  the 
estimation  of  the  energy  spectrum  in  the  form  defined  here. 

More  important,  thoughts  the  illumination  shed  on  the  problem 
of  formulating  signal  processing  systems  intended  for  use  v/ith  non¬ 
stationary  backgrounds  or  in  situations,  such  as  explosive  echo 
ranging,  where  the  waveform  to  be  detected  is  itself  a  random  trans¬ 
ient.  We  have  seen  that  under  certain  conditions  of  envelope  duration 
and  noise  characteristic  that  it  is  possible  for  the  spectra  and 
autocorrelations  of  the  individual  transients  to  be  rather  different 
at  the  same  value  of  an  argument,  even  though  the  expectations  are 
the  same  from  sample  to  sample.  This  implies  immediately  that  it  may 
not  be  advisable  to  tailor  the  characteristics  of  a  monopulse  process¬ 
ing  system  too  closely  to  the  mean  values  of  autocorrelation  and 
spectrum.  A  filter  painstakingly  devised  to  reproduce  or  complement 
the  mean  spectral  density  may  well  do  serious  violence  to  the 
individual  transients  just  because  they  very  well  may  have  spectra 
which  differ  significantly  from  the  mean.  The  same  considerations 
apply  to  correlation  processing.  It  is  hoped  that  the  present  theory 
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contributes  at  least  some  understanding  of  this  problem  and  aids  in 
quantifying  the  latitude  that  must  be  allowed  for  accommodating  pulse 
to  pulse  differences.  Admittedly,  a  rigorous  treatment  of  this 
question  must  await  the  extensions  described  above,  particularly  the 
derivation  of  joint  and  conditional  probabilities  and  a  suitable 
linear  system  theory.  These  first  considerations,  though,  should  at 
least  alert  the  researcher  to  the  existence  of  the  problems  involved. 
It  is  especially  hoped  that  more  care  will  be  taken  in  the  measure¬ 
ment.  and' description  of  such< phenomena  as  sonar  reverberation  in 
connection  with  the  study  of  monopulse  detection  systems.  There  is 
little  evidence  that  these  areas  have  been  approached  in  the  past 
with  the  rigor  they  deserve. 


EDWARD  C.  WHITMAN 

Magnetics  and  Electrical  Division 
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APPENDIX  A 

ON  BOUNDING  THE  INTEGRAL  j/<Ppp(u,o)  V^n+r)  ^n(n-r)du 


Since  the  integrand  is  an  even  function,  we  can  consider  the 
infinite  integral 


i^(r)  =  /Vu,o)^-r)  ^nn(u-r)du 

-00 


and  write  it  as  follows: 


A-l 


ix(T)  =/V^p(u,o)  ^nn(u+r)  ^pp(u,7)  VnnCu-^du 
—  ©0. 

Now  using  the  Schwarz  inequality  as  in  reference  (e) 


v°° 

hH2- 

/^ppCu,0)  ^(u+r)du 

^  -  0© 

J  ‘fppCujo)  ^nnCu-rJdu 

m 

A -2 


A-3 


Since  VU’°?  is  a  real  even  function,  the  factors  on  the  fight  are 
equal  and  positive.  Therefore, 


ii(r) 


< 


( u ,  6 )  'V,n^(u-T)du 


h-k 


Consider  now  the  integral  on  the  right  of  A-4  arid  denote  it  i?(T) . 
Thus, 


il(r)  ~  i2(r) 

Since  is  a  real 

therefore 


=  ^(u-l^du 

—  OO 

even  function,  l^n^Cu-T)  r  ‘V/J:(r-u) 


A-5 


and 
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/%  Od 

i2(^)  *  pp(n»0)  ^(r-u)QU  =  <fpp(u,0)  9  ^(u)  A-6 

—  PO 

where  8  denotes  the  operation  of  convolution.  Row,  let  us  enter 
the  frequency  domain  where  the  transforms  of  i2(r),  ?pp(uj°)i  and 

V2(u)  are  respectively  P(w),  and  Q(<o).  In  reference  (a) 

it  is  pointed  out  that  the  Fourier  transform  of  a  realizable  auto¬ 
correlation  function  must  always  be  positive.  Thus  P-.(»)  is  always 
positive,  and  since  Q(«)  is  the  transform  of  the  square  of  a  real¬ 
izable  autocorrelation  function,  it  must  be  the  convolution  of  a 
positive  transform  with  itself.  Thus  Q(«)  is  also  everywhere  greater 
than  zero.  Writing  (A-6)  in  the  transform  domain  yields 

I2(w)  =  P(«)  Q(«)  A -7 

and  re-transforming  to  the  time  domain  gives 


/oo 

P(«)  Q(»)e^  “T  du 


i2(r)|s/TpMlQ(«)U«  , 


but  since  by  A-6,  i2(*)  must  always  be  positive  and  since  P(»)  and 
Q(\«)  are  greater  than  zero,  this  becomes 


OO 

i2('r)<  f  p<«)  'Kw)d«  *  i2(0) 


A -10 


By  A-5, 


r)  <  i2(r)  <  i2(0) 


A-ll 


and  finally 


Wu-^d«</~pp<u>0)Tm<u>  dU 


A-12 
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when  ^>pp(]a,0  )  is  always  positive.  Under  this  restriction,  it  has 

been  sho’-'n  that  the  maximum  of  the  integral  of  interest  occurs  when 

T  s  0. 


A -3 
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APPENDIX  B 

AN  ALTERNATIVE  APPROACH  TO  THE  NOISE  BURST  SPECTRA  USING  THE 

WIENER -KHIKTC H I NE  RELAT I ON 


As  applied  to  transient  signals,  the  Wiener-Khintchirie  relation 
states  that  the  autocorrelation  function  of  a  signal  (as  defined  in 
equatiom  (4-).  of  the  main  body)  and  its  energy  spectrum  constitute  a 
Fourier  trans/form  pair: 


^ss^ 


■p 


(w)<r 


B-la 


♦Moo  *  rb*  J  ^ss^)e“jMrdr  B-ib 

-oo 

Thus,  it  appears  that  one  can  compute  the  energy  density  spectrum 
from 'the  second  of  -these  equations,.  From  our  previous  considerations, 
however,  we  know  that  <£„,.(«):  is  a  random  variable  for  all  u>  and 

that  we  thus  must  be  content  with<  computing  means  and  variances. 


e[*ss(»)]  * 


By  equation  (23)  of  the  main  body, 

eBTssW]  = 

and  thus  ^ 

E^ss^D  =  -j? j  *nr,W  B-3 

-  OO 

precisely  the  same  expression  found  in  the  main  body  from  the  Fourier 
integral  approach  (equation  (64)). 


E-l 
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Turning  now  to  a  calculation  of  the  variance  of  4>ss(to)  from 

this  standpoint,  we  must  first  compute  the  mean  square  value  of  the 
spectrum  for  each  a>  . 


2  = 


s.<*>*,“r  d 


’  ^fj' 


ss(r)  (T+tr)dr  d‘ 


-CA  - Oo 


£..K]  =  ff  E  [♦«(*>  ^S^J 


-Mr+<r) 

e  d  r  d  *r 


— oo  —  o* 


Mow  by  definition, 


./: 


^'ss^J  =  /  e(t )e(t+r)n(t )n(t+T)dt  /  e(x)e(x+<r)n(x)n(x+<r)dx 


#• 


t)e(x)e(t+r)e(x+<r)nCt)n(x)n(t+r)n(x+<r)dt  dx 


woo  -C» 


Oo  OOl 

S&..W  ‘fcss^  zff )e(x)f\(t+r)e(x+<r)  E[n(t)n(x)n(t+r)n(x+o^]dt 


and  combining  (B-5)  and  (B-7),  there  emerges  that 


09  an  oo  orz 


=  ■jj^r ffff c(t)e(x)e(t+r)e(x+«r)  E[n(t)n(x)n(t+r)n(x+tf"T] 


_C*»  o  o 


3*.  <53-  ^-8 
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Now  with  the  substitution  of  variables  us  t+T,  v  ■  x  +<r  , 


(t)e(x)e(u)e(v )  E[n(t)n(x)n(u)n(V/7] 


r*a  (u-t+v-x)dt  dx  ^  dv  B-9 


which  is  equivalent  to  the  integral 


OO  <JW  (M 

=  tk  f/ff  e(t)e(x)e(u)e(v)  E[n(t)n(x)n(u)n(v)] 

e-jo>  (t-x+u-v)  dt  dx  du  dv  B-10 


o  o  0  » 


This  expression  is  identical  with  that  found  as  equation  (73)  in 
approaching  the  variance  from  Fourier  integral  considerations.  If 
both  the  mean  and  mean  square  are  identical  for  the  two  methods j  the 
variance  must  be  also  and  hence  it  appears  that  the  WienemKhintchine . 
relation  is  directly  applicable  to  the  expectations  of  the  statistical 
quantities  associated  with  random  transients. 
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